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Development  of  a  finite  volume  code  for  spatially  evolving  simulations  of  flow  past  wings: 

We  have  developed  a  three-dimensional  finite  volume  code  for  performing  simulations  of  spatially 
evolving  flows.  The  finite  volume  technique  is  preferred  due  to  its  flexibility  for  implementing  a  variety 
of  realistic  and  physically  relevant  boundary  conditions  e.g.  unlike  spectral  codes  this  method  is  not 
constrained  to  periodic  boundary  conditions  only  and  therefore  allows  simulations  of  flows  past  wings 
with  finite  spans.  The  algorithm  is  formulated  with  equations  written  in  generalized  coordinates  which 
makes  this  code  useful  for  simulating  a  wide  variety  of  flows  including  spatially  evolving  jets,  flows  past 

bluff  bodies  of  different  geometries,  etc.  Also,  versatile  grid-generation  is  useful  in  compressible  flows 
for  local  grid  refinement  to  resolve  shocks. 

Grid  generation  is  performed  m  a  preprocessor  routine  which  generates  a  body-fitted  ,  three  dimensional , 
orthogonal  grid  system.  The  code  computes  the  Cartesian  coordinates  of  grid  points  in  the  transformed 
system  and  their  corresponding  scaling  factors  to  be  accessed  by  the  main  program.  This  allows 
versatihty  in  the  specification  of  the  inner  boundary  which  could  correspond  to  an  airfoil,  a  cylinder, 
etc..  The  Cartesian  system  is  used  in  the  computational  domain,  while  the  generalized  system  is  the 
physical  domain.  Once  the  grid  is  generated,  the  Jacobian  of  the  transformation  (from  the  physical 
domain  to  the  computational  domain)  is  evaluated  at  eveiy  point.  From  the  value  of  the  Jacobian,  the 
scaling  factors  are  calculated.  A  noniterative  method  of  grid  generation  using  parabolic  differential 
equations  is  preferred  .  The  advantages  of  using  a  set  of  parabolic  equations  are: 

•  Parabolic  equations  are  initial  value  problems,  so  grids  are  generated  by  a  marching 
algorithm  like  the  hyperbolic  grid  generation  method. 

•  The  parabolic  partial  differential  equations  have  most  properties  of  the  elliptic  equations, 

particularly,  the  diffusion  effect  which  smoothes  out  any  singularity  of  the  inner  boundary 
condition. 

•  By  employing  the  marching  algorithm,  the  computational  time  required  is  only  a  very  small 
fraction  of  that  for  the  elliptic  grid  generation. 


Tl.=  govsnimg  equations  are  formulaied  using  finite  volume  discreUsation ,  with  velociUes  and  ptessure 
as  the  pnmary  variables.  The  velocity  and  pressure  points  are  distributed  over  a  staggered  grid.  The 
algorilhm  involves  a  ftactlonal  step  nwUtod  ( predictor  step  and  conecrer  step )  to  stabilize  the  advective 
tenns.  The  oKfitod  involves  an  explicit  evaluation  of  dte  velocities  in  the  predictor  and  corrector  stitges 
and  an  implicit  evaluation  of  the  auxiliary  pressure.  The  code  is  second-order  accurate  in  space  and  time 
(the  code  has  been  implemented  in  a  manner  which  allows  easy  extension  to  a  compaet-siorage  third 
order  Runge-Kuita  scheme);  a  scheme  similar  to  the  SMAC  method  is  used  for  solving  the  Navier- 
Stokes  equations.  We  have  implemented  this  algorithm  and  validated  it.  Validation  cases  include  iwo- 
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dimensional  flow  in  a  square  lid-driven  cavity,  using  a  uniformly  spaced  Cartesian  grid.  Streamlines, 
vorticity  contours,  center-line  velocity,  peak  vorticity,  etc.  obtained  from  this  simulation  were  found  to  be 
in  good  agreement  with  the  results  in  the  literature.  In  Figs.  6,7  we  show  vorticity  contours  in  a  lid-driven 

flow  in  a  square  cavity  for  two  different  Reynolds  numbers,  these  compare  very  well  with  documented 
results  in  the  literature, 


(In  the  following,  .superscripts  denote  the  relevant  publication  from  the  list  at  the  end  of  the 
report) 

A  model  of  turbulent  vortex  breakdown^: 

We  have  developed  a  theory  of  turbulent  vortex  breakdown  to  predict  hysteretic  jump  transitions  in 
recirculatory  2one.s  over  delta  wings  and  to  explain  recent  experimental  observations  of  conical  vortex 
breakdown.  Our  approach  involves  a)  Modeling  swirling  jets  by  conically  similar  solutions  of  the  Navier- 
Stokes  equations  (Fig.  1  shows  the  destroyed  vortex  in  spherical  coordinates  (r,0,())),  curves  3.4  show 
typical  .streamlines  of  the  meridional  motion  in  the  two  flow  cells  separated  by  a  conical  surface),  b) 
modeling  turbulence  by  an  eddy  viscosity  which  is  represented  by  a  step  function  of  the  polar  angle  and 
c)  matching  solutions  for  the  laminar  and  turbulent  flow  regions.  The  model  shows  the  formation  of  a 
conical  turbulent  wake  which  develops  without  a  recirculatory  zone.  Fig.  2  shows  the  results  of  the  these 

calculations  where  the  velocity  at  the  axis,  v,  is  plotted  as  a  function  of  the  long  parameter,  M.  This 
parameter  is  defined  as 

M  =  21tJo/Re.^ 

where,  Jo  is  a  measure  of  the  flow  force  acting  on  the  plane  normal  to  the  axis  and  Re,  is  the  Reynolds 
number  of  the  vortex  based  on  its  circulation.  The  plot  shows  the  velocity  on  the  axis  for  both  turbulent 
(T)  and  laminar  (L)  flows.  Positive  values  of  correspond  to  the  consolidated  jet  and  negative  values  to 
the  two-cell  flow.  One  can  see  that  the  curves  have  folds  at  smaller  values  of  M=Mf  than  M=Ms  for  flow 
separation  (v„=0).  The  results  obtained  for  the  laminar  flow  coincide  with  those  obtained  by  Long.  Our 
results  also  provide  upper  estimates  for  the  angles  of  the  conical  wake  region  where  the  turbulence  is 
strong  and  for  the  vortex  core  where  the  swirl  is  of  the  solid  body  type.  The  prediction  of  the  conical 
wake  and  the  estimates  for  its  angle  also  agree  with  earlier  studies. 

Vortex  breakdown  control  using  vortex  generators^: 

To  develop  effecUve  means  to  control  vortex  breakdown  in  wins-dp  and  leading  edge  veniees.  we  have 
studted  swirling  flows  in  a  conflned  cylinder.  Such  a  confined  flow  is  relevan,  to  the  study  of  this 
probien.  as  tt  ts  free  from  ambient  disturbances.  The  basic  flow  is  driven  by  a  rotating  disc  and  is 
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characterized  by  a  Reynolds  number,  Reo.  A  thin  vortex  core,  which  develops  at  the  axis  at  large  Rcd. 
undergoes  breakdown.  The  core  is  the  most  sensitive  flow  region  and  is  forced  by  an  independently 
rotating  rod  (with  Reynolds  number,  Rcr  -  Rep  /  100)  which  has  a  diameter  comparable  to  that  of  the 
core.  Flow  visualization  revealed  that  even  weak  forcing  by  the  rod  can  strongly  influence  the  flow  (Fig. 
3  ).  With  increasing  rod  co-rotational  speed,  the  size  of  the  breakdown  bubbles’  decreases  and  eventually 
a  conical  wake  without  any  recirculating  regions  is  formed.  In  contrast,  when  the  rod  is  counter-rotated, 
the  number  and  size  of  the  bubbles  increases  (Fig.  4).  We  have  developed  a  theoretical  model  to  explain 
the  physics  of  these  flow  pattern  transformations.  We  have  found  that  the  flow  with  a  co-rotating  rod  can 
be  modeled  using  conically  similar  solutions  of  the  Navier-Stokes  equations.  Fig.  5  shows  the  modeled 
shape  of  the  wake  when  vortex  breakdown  is  suppressed  which  closely  corresponds  to  the  wake  surface 
observed  via  flow  visualization. 

Inviscid  Vortex  Breakdown^: 

We  demonstrated  that  axisymmetric  vortex  breakdown  in  a  steady  flow  in  a  semi-infinite  circular  pipe  is 
a  reversible  process  which  appears  due  to  solution  nonuniqueness  in  some  range  of  inflow  parameters, 
when  the  entire  flow  jumps  to  another  metastable  steady  stale  with  the  same  boundary  conditions.  All 
such  solutions  are  smooth  along  the  pipe  length  and  have  the  same  mechanical  energy,  but  the  flow 
forces  in  the  different  flow  states  are  different.  Vortex  Breakdown  must  appear  by  continuous  change  in 
flow  parameters  (usually  the  swirl  number)  when  locally  the  solution  fails  to  exist  because  of  fold  or 
similar  catastrophe  which  can  be  related  to  transcritical  bifurcation  of  a  degenerate  state  (characterized 
by  the  existence  of  trivial  solutions),  but  the  spontaneous  jumps  between  different  solutions  (not  on  a 
fold)  also  are  considered  as  a  finite  perturbation  instability.  Among  the  possible  ways  to  eliminate 
degeneracy,  we  consider  the  injection  of  azimuthal  vorticity  into  the  vortex  core  at  the  pipe  entrance.  We 
show  that  the  flow  pattern  inside  the  steady  near-axis  separation  zone  (usually  associated  with  vortex 
breakdown)  depends  on  the  flow  history.  Therefore,  our  approach  to  vortex  breakdown  takes  into  account 
the  flow  pattern  inside  the  separation  zone.  We  consider  two  approaches:  the  traditional  analytic 
continuation  method  (leading  to  a  recirculation  zone)  and  the  stagnant  separation  zone  model  (without 
velocity  jump  on  its  boundary).  We  reveal  serious  defects  of  the  analytic  continuation  method:  infinite 
negative  velocity  at  a  finite  inflow  swirl  and  global  nonexistence  of  solutions.  Analytic  continuation  can 
only  lead  to  jump  transition  from  a  stable  state  to  an  unstable  state  or  internal  separation,  but  not  vortex 
breakdown.  The  new  stagnation  zone  model  is  superior  in  that  solutions  always  exist  and,  for  large 
enough  inflow  swirl,  exhibit  nonuniquenes  and  folds  by  smooth  variations  of  the  flow  parameters,  thus 
predicting  the  experimentally  observed  hysteretic  jump  transitions  in  vortex  breakdown.  This  model  has 
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revealed,  m  addition,  some  new  phenomena,  such  as  vortex  contraction  and  a  finite  distance’ effect:  with 
increasing  swirl  a  stagnation  zone  of  nonzero  size  suddenly  appears  far  downstream  and.  consequently, 
the  stagnation  zone  inside  the  pipe  abruptly  arises  at  a  finite  distance  from  the  inflow. 

Vortex  Breakdown  as  inviscid  separation  in  steady  flows^: 

We  swdied  vortiej  channel  flows  with  different  inflow  vorticity  distributions  and  showed  that  an  exactly 
stagnant  separation  none  is  formed  in  the  linht  that  the  inflow  vorticity  becomes  discontinuous.  It  is 
shown  that  the  analytical  continuation  model  applied  to  taotarional  channel  flow  with  a  slender  vortical 
cote  (the  planar  analog  of  the  RanKne  vortex)  yields  Vorticity  bteakdown'-a  phenomenon  similar  to 
vortex  breakdown  in  swirling  flows.  This  study  is  related  to  the  classical  problem  of  the  indeterminacy  of 
steady,  inviscid,  incompressible  vortical  planar  flows  due  to  presence  of  separation.  We  have  proved  that 
any  inviscid  steady  flow  is  'totally  nonunique'  in  the  sense  tha.  there  exists  a  continuum  of  possible 
solutions  mvolving  tecirculadon  zones  which  satisfy  the  same  boundary  conditions.  For  sufficiently  large 
inflow  vorUcity,  the  formation  of  separation  zones  is  inevitable.  Various  additional  phenomenological 
conditions  to  achieve  uniqueness  and  to  determine  the  flow  in  the  separation  zone  were  evaluated.  We 
revealed  a  new  'measure  paradox':  a  stagnation  zone  of  finite  area  arises  from  one  streamline  by  changing 
vorticity  at  only  one  inflow  point.  We  focus  on  this  new  steady  stagnant  zone  with  continuous  velocity 
distributiou  everywhere.  The  existeuce  and  uniqueness  theorems  for  a  flow  with  sragnaUon  zone  were 
proved.  We  showed  that  the  kinetic  energy  of  the  entire  flow  domain  is  minimized  when  the  separation 
zone  is  stagnant  and  is  completely  contained  within  any  possible  tecirculadon  zone.  Nonuniqueness  and 
hysteresis-related  multiple  solutions  are  inherent  to  those  solutions  with  recirculation  zones  where  very 
large  (physically  unrealisUc)  velocity  in  the  reversed  flow  appears.  Being  free  from  these  drawbacks,  the 
stagnation  zone  model  appears  to  be  superior.  A  numerical  investigation  of  a  time-evolving  voriical  flow 
m  a  square  domain  confirmed  the  existence  of  a  critical  value  of  inflow  vorticity,  above  which  the  flow  is 
unsteady,  depending  on  the  initial  condition  within  the  domain. 

Hysteresis  in  swirling  jets^: 

Vortex  breakdown  can  be  viewed  as  a  hystereHc  transformation  in  swirling  jet.  We  consider  viscous 
flows  wrrh  steady  rotationally  symmetric  motion  above  an  infinite  conical  stream  surface  with  a  half¬ 
angle.  e.  The  flows  analyzed  ate  generalizations  of  Long's  voncx.  They  correspond  to  conically  similar 
solutions  of  the  Navier-Stofces  equadons  and  are  characterized  by  circuladon,  T,  given  at  the  surface  and 
axial  force  I, .  Asymptodc  analysis  and  numerical  computadons  show  that  four  ( for  %  S  90>)  or  five  (6,  > 
90")  soludons  exist  in  some  range  of  T,  and  I,.  The  solution  branches  form  hysteresis  loops  which  ^ 
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related  to  jump  transitions  between  various  flow  regimes.  Four  kinds  of  jumps  were  found;  (i)  vortex 
breakdown  which  transforms  a  near-axis  jet  into  a  two-cell  flow  with  reverse  flow  near  the  axis  and  an 
annular  jet  fanning  out  along  conical  surface  9  :=  0,  <  0,;  (ii)  vortex  consolidation  causing  a  reversal  of 
(i);  (iii)  jump  flow  separation  from  surface  0  =  0,;  and  (iv)  jump  attachment  of  the  swirling  jet  to  the 
surface.  As  T,  and  /  or  J,  decrease,  the  hysteresis  loops  disappear  through  a  cusp  catastrophe.  In 
particular,  vortex  breakdown  can  be  interpreted  as  a  fold  catastrophe.  Two  new  and  striking  effects  were 
found;  (i)  there  is  a  pressure  peak  of  Offc^)  inside  the  annular  swirling  jet;  and  (it)  a  consolidated 
swirling  jet  forms  with  a  reversed  (’anti-rocket’)  flow  force. 
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Fig  1 :  Schematic  of  the  problem 
and  the  coordinates 
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Fig  2:  Fold  catastrophe  in  the  laminar 
(L)  and  turbulent  (T)  flows 


Fig.  3;  Flow  visualization  pictures  and  schematics  of  streaklines  for 
co-rotating  disc  and  central  rod.  Rejj=2720.  (a)  Rer^O;  (b)  Rer=21;  (c) 

Re|.=29.  Arrows  in  the  schematics  indicate  flow  direction. 
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Fig  4:  Flow  visualization  pictures  showing  time  evolution  of  multiple 
breakdown  bubbles  for  the  case  of  counter-rotating  disc  and  central  rod.  Recj- 

2720.  Rer=14.5. 


Fig.  5:  Analytical  solution  for  the  shape  of  the  conical  wake  observed  in  the 
case  of  co-rotating  disc  and  rod. 


Fig.  6:  Vorticity  contours  in  a  lid-driven  cavity 
flow  at  Re=100:  t=1.5 


Fig.  7:  Vorticity  contours  in  a  lid-driven  cavity 
flow  at  Re=400;  t=9 
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